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Abstract. When subjected to a linearly polarized terahertz pulse, a mesoscopic ring 
endowed with spin-orbit interaction (SOI) of the Rashba-Dresselhaus type exhibits 
non-uniform azimutlial charge and spin distributions. Both types of SOI couplings 
are considered linear in the electron momentum. Our results are obtained within a 
formalism based on the equation of motion satisfied by the density operator which 
is solved numerically for different values of the angle 4>, the angle determining the 
polarization direction of the laser pulse. Solutions thus obtained are later employed in 
determining the time-dependent charge and spin currents, whose values are calculated 
in the stationary limit. Both these currents exhibit an oscillatory behavior complicated 
in the case of the spin current by a beating pattern. We explain this occurrence on 
account of the two spin-orbit interactions which force the electron spin to oscillate 
between the two spin quantization axes corresponding to Rashba and Dresselhaus 
interactions. The oscillation frequencies are explained using the single particle 
spectrum. 
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1. Introduction 

The controlled induction of spin and charge currents in mesoscopic rings represents a 
longstanding interest of experimentalists and theorists alike. By taking advantage of 
the right-left asymmetry of the electron states, various methods have been proposed 
over the years to this end [H El [3]. In particular, by applying an ultrashort, terahertz 
frequency laser pulse with spatial asymmetry, the charge and spin excitations are 
realized in a non-adiabatic regime, since their excitation occurs on a time scale that 
is much shorter than the electron relaxation lifetime [21 HI [5]. More recently, these 
ideas have been reconsidered in the case of rings made out of zinc-blend semiconductor 
structures that exhibit spin-orbit interaction (SOI). The linear coupling between the 
electron momentum and spin appears as a result of broken inversion symmetry that 
results from spatial confinement (Rashba - R) or from bulk properties (Dresselhaus 
- D) [7]. Under these general circumstances, it was shown that the spin and charge 
flow simultaneously [HI |9l HO] and special configurations, such as hybrid structures [H] 
or a number constraint (odd number of particles) [12] are needed to isolate the spin 
contribution. Recently, we demonstrated that by applying a two-component terahertz 
frequency pulse, for certain values of the phase difference angle <^ between the two 
components, a pure spin current is non-adiabatically induced in the ring. Physically, 
this situation is a consequence of the interplay between SOI which rotates the electron 
spin around the ring and the asymmetry of the external perturbation embodied by 0. 

The present paper discusses a follow-up to this analysis by studying the effect of 
a terahertz laser pulse on a ring endowed with both Rashba and Dresselhaus spin- 
orbit couplings. The coexistence of these couplings introduces precession effects of 
the electron spin along different directions, a phenomenon bound to affect the spin 
and charge distributions around the ring. Our formalism is based on obtaining the 
stationary solutions of the equation of motion satisfied by the particle-density operator 
which are later employed in estimating the spin and charge currents as as statistical 
averages. The three independent parameters of the problem are the phase of the laser 
excitation, 0, and the coupling strengths of the Rashba and Dresselhaus interactions. 
Our main result is that, for comparable R and D intensities, the spin current, calculated 
along the z direction, showcases an oscillatory pattern that results from the nutation of 
the electron spin between the two natural directions associated with the two spin-orbit 
interactions. 

2. The Equilibrium Hamiltonian 

We consider a one- dimensional (ID) quantum ring of radius tq in the presence of both 
Rashba and Dresselhaus spin-orbit interactions. The ring is discretized in sites with 
a lattice constant a = 27rro/N. The polar coordinates are chosen such that the site with 
index n has an azimuthal angle 6n = 2mi/N . We denote by c^^ and Cn^a the creation 
and annihilation operators for an electron with spin a on the site n = 1, A^. 
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Figure 1. The energy spectrum of a ring with SOI. The Rashba energy Va = 0.1 and 
the Dresselhaus energy Vg = 0.05. On the horizontal axis is the rank of the energy 
value within the spectrum of the Hamiltonian. All states arc spin degenerate. The 
energy unit is 1/ = fi,^/2m*a^. At the initial time only the first 6 energy levels are 
occupied. 

In the following considerations we adopt the discrete representation of the 
Hamiltonian -faring = Hq + Vr + Vd , that describes the non- interacting electron dynamics 
around the ring in the presence of the Rashba (R) and Dresselhaus (D) SOI couphngs. 
Its components are: 

Ho = 2VY^ C+^C^^^ -VY^ C+aC„+l,^ -^Yl Cn,aCn-l,<x, (l) 
n,a n,a n,<T 

Vr = - iV^ Yl ^r{en,n+ir''''4.C^+l^' + h.C. , (2) 
n,a,a' 

Vd= - iVp Y ^e*(^n,n+l)"''^'ctc„+i<,, + h.C. . (3) 
n,a,(T' 

Three different energy scales are introduced: the hopping energy V = h?/2m*a'^ , 
Rashba energy Va = a /2a and Dresselhaus energy Vg = P/2a, with a and /3 being 
SOI parameters and m* the effective mass of the electron in the host semiconductor 
material. It is customary to introduce the radial and azimuthal spin matrices cxr and ag 

ar{0) = a.j. cos 9 + ay sin 9 , ag{9) = —a^,; sin 9 + ay cos 9 . (4) 

In the presence of only one type of SOI, either Rashba or Dresselhaus, the eigenstates 
of the Hamiltonian can be obtained analytically [9l [101 Ej- In both cases, the spin 
quantization axis is tilted relatively to the z axis by angles Qa and Qjs respectively, 
given by 
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where A6' = 2tx/N. Consequently, if either = or = 0, there exist a 
corresponding natural direction of the spin, given by = cosfi^e^ — sinfio-e,, or 
= cosf2^e^ + sinfi/je.e, respectively, along which the projection of the electron spin 
is conserved and equal to ±/i/2. The directions and are called the principal spin 
axes, or Rashba and Dresselhaus spin quantization axes, respectively. The difference 
between the Rashba and the Dresselhaus SOIs is reflected by the orientation of the spin 
precession relatively to the orbit: in the Rashba case, the precession and the orbital 
angular speed have the same sign, whereas in the Dresselhaus case they are opposite 

[131 do]. 

When both interactions are present the diagonalization of the Hamiltonian can be 
performed only by numerical methods. The energy spectrum for = 20 sites and 
with Va = O.IV and Vg = 0.05V is shown in Fig. [H Each state is double (spin) 
degenerated. For a ring radius ro = 50 nm our SOI strengths become comparable to 
those in experimental materials [13]. One obtains a ~ 21, (] ^ 10 meVnm for InAs (i. 
e. using m* = 0.023me), ot a ^ 34, /3 f» 17 meVnm for InSb (i. e. using m* = 0.014me). 
Our time unit is t = h/V, corresponding to ^ 0.1 ps for InAs and ~ 0.06 ps for InSb. 
We will consider rig = 6 (noninteracting) electrons in the ring. The occupied states have 
energies E12 = —0.01184, £'3 4 = 0.07311, -Es^g = 0.1003. The next three energy values 
are Ej^s = 0.3484,^9,10 = 0.3967, and En,i2 = 0.3967 (units of V"). The even and odd 
labels correspond to positive and negative spin chiralities, respectively. 

The combination of the two types of SOI breaks the circular symmetry of the 
electron distribution and creates a charge deformation as shown in Fig. O This effect has 
been discussed in the recent literature both for one-dimensional [10] and two-dimensional 
[T4] quantum rings. In the one-dimensional ring model the charge density has two 
minima along the direction y = x and two maxima along y = —x. The effect can be 
explained by an intrinsic periodic potential that develops around the ring |10j . 

U{e) = 2VaVfism{2e). (6) 
3. Time dependent perturbation 

At time t = the quantum ring is excited by a short terahertz laser pulse linearly 
polarized along a direction making the angle with the y axis. The pulse has a lifetime 
of l/F and an angular frequency u and can be represented by the time dependent 
Hamiltonian 

Hpn\se{t) = Ae~^* sin(a;t) cos(6' + </>) . (7) 

The time evolution of the system is given by the Liouville equation, 

ihp{t) = [H + Hp^Ut), p{t)] , (8) 

where p{t) is the time dependent statistical operator. The time evolution of the electron 
states, assumed to be initially in the minimum energy configuration, is described by the 
numerical solutions of the Liouville equation obtained with the Crank-Nicolson finite 
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Figure 2. The electron density around the ring. Due to the presence of both SOIs 
the distribution of electrons is nonuniform. Minima occur along the direction y ~ x, 
or at angles 6 = tt/A and 9 = 57r/4, and maxima along y = ~x or at angles 9 = 37r/4 
and 9 = T-k/A. The blue dots indicate the = 20 sites on the ring. 



difference method [3] with small time steps (5t ^ F ^. The expectation value of any 
observable O is then calculated as (O) = Tr [p{t)0]. 

4. Numerical results 

In the following example the pulse parameters are A = 2.3V, u = 0.3V/ h, and 
r = cj, generating a pulse lifetime of 3.3h/V. The time evolution of the charge current 
Ic{t) = ev{t) for = is shown in panel (a) of Fig.[3l After the pulse vanishes, the charge 
current exhibits a large period oscillatory behavior. At high frequencies, superimposed 
secondary oscillations are noticeable. These findings can be explained on account of 
single electron excitations. The large period is correlated with transitions between the 
second and the third states on the energy scale, {3,4} and {5,6} respectively, that 
generate 2tt/{Eq — E4) ^ 231h/V. Transitions between the states {3,4} and {1,2}, 
occur with period Tih/V, i. e. one third of the large period, have less weight and cannot 
be distinguished from the large oscillations. Excitations over larger energy intervals, 
like 2-71/ {Es — Eq) ^ 23h/V, or 2n/{Eio — E2) ~ Ibh/V, etc. are responsible for the 
finer oscillations, with periods varying from bOh/V down to I'bh/V or less. 

Next, in Fig. [3](b) we illustrate the time evolution of the spin current, defined by 
-^Sz — 1 {^z^d + vgaz) , i. e. polarized along the z direction, for = 0. This instance 
indicates more clearly that the time evolution of the system is not periodic even long after 
the external pulse has vanished. A striking beating pattern is observed, on a background 
of fast oscillations, faster and much more pronounced than for the charge current. To 
explain the fast oscillations we refer to Fig. HI which shows the time dependence of the 
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Figure 3. Charge and spin currents in the presence of both Rashba and Dresselhaus 
SOI. Va = O.IV and Vg = 0.051^. The upper panels show the resuhs for the relative 
angle = 0, and the lower panels for </> = 7r/5. The oscillations are explained by the 
Bohr frequencies of the energy spectrum, see text. The green lines in the upper panels 
show the results obtained in [5] with only one type of SOI interaction. The time is in 
units of h/V 



spin projection on the z axis within a time interval of 200 units. A series of double peaks 
is noticeable created by the two SOI couplings that generate spin torque in opposite 
directions. The result is a superposition of Rashba and Dresselhaus precessions, with a 
certain relative phase shift. Consequently, the fast oscillations seen in the charge current 
are doubled for 5*^ and for the the spin current. The average time between consecutive 
maxima in Fig. H] is about 7.4 h/V which is about half of the time corresponding to the 
fast oscillations of the charge current within the same time interval. 

We find, therefore, that the spin current is much more sensitive to the single electron 
transitions between states with distant energies, implicating high frequencies, than the 
charge current. The beating pattern is a result of that. Since the SOI parameters are 
different, Va 7^ Vg, the principal spin axes and e/3 have different angles with the z 
direction, fi^ 7^ fi^, and the electron spin executes rapid oscillations between them. In 
physical terms the orbital motion is accompanied by spin nutation. 

At equal SOI parameters, Va = Vp, the principal directions do not coincide, but 
make the same angle with the z axis, Via = ^13- The spin precessions cancel each 
other and the corresponding spin current vanishes. However the charge current remains 
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Figure 4. The spin projection Sz witliin a time subinterval. The fluctuations look 
like a series of split peaks. The time is in units of h/V 




Figure 5. Amplitude of the charge current vs relative angle 



qualitatively as in Fig. |3l If only one type of SOI is present, the spin current is not zero, 
but constant, whereas the charge current vanishes [5], as shown by the green lines in 
Fig. [3] (a) and (b). 

In the lower panels of Figs. [3l (c) and (d), we show the results for the charge 
and spin currents calculated at (p = tt/S. When compared with the previous case, the 
amplitude of the slow oscillations of the charge current decreases, while that of the spin 
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current increases. Fast oscillations are still pronounced in the spin current. Moreover, 
the spin current shows similar oscillations as the charge, with a phase shift of 7r/2, 
which dominate now the nutation effects. This effect can be explained on account of the 
anisotropy of the charge distribution shown in Fig. |2l The symmetry axis of the ring 
corresponds to ^ = 7r/4, and therefore if the pulse is oriented along this direction no 
charge current is induced. This is seen further in Fig. [5] where we show the amplitude 
of the charge current vs. the angle 0, which vanishes for = 7r/4. In this case a spin 
current of maximum amplitude is induced due to the opposite chiralities of the states. 
The slow oscillations overlap with the fast fluctuations described earlier. At = 0, the 
charge current has a maximum amplitude, whereas that of the spin current is minimal 
and reduced to the fast fluctuations related to spin nutation. The angle </> = vr/S is used 
in Fig. [3]^c) and (d) as the closest value to 7r/4 which can be obtained on our discretized 
ring with 20 sites. 

5. Conclusion 

The simultaneous effect of Rashba and Dresselhaus spin-orbit interactions upon spin 
and charge currents are estimated numerically in a ring geometry. A linearly polarized 
terahertz laser pulse is used for a non-adiabatic generation of the currents. In the 
presence of the dual spin-orbit coupling, the charge and spin currents exhibit a quasi- 
periodic behavior whose amplitude depends on the relative angle between the pulse 
and the charge deformation created by the SOI effects. The spin current reflects a 
nutation of the electron spin between the Rashba and Dresselhaus spin quantization 
axes, superimposed as a fast beating pattern on the main oscillation. 
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